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iCJ ■ Abstract 

The effective action for the Polyakov loop serving as an order parameter for 
deconfinement is obtained in one-loop approximation to second order in a derivative 
expansion. The calculation is performed in d > 4 dimensions, mostly referring to 
the gauge group SU(2). The resulting effective action is only capable of describing a 
deconfinement phase transition for d > d CI ~ 7.42. Since, particularly in d = 4, the 
system is strongly governed by infrared effects, it is demonstrated that an additional 
infrared scale such as an effective gluon mass can change the physical properties of 
the system drastically, leading to a model with a deconfinement phase transition. 

• rH 

1 Introduction 

As a prelude to a truly nonperturbative evaluation of the effective action of Yang-Mills 
theory, the one-loop effective action with all-order couplings to a specific background may 
provide a first glance at the up to now unknown ground state of the theory. Since the 
problem of confinement is supposed to be intimately related to the quest for the ground 
state, it is elucidating to investigate the response of several "confining vacuum candidates" 
on quantum fluctuations - even in a perturbative approximation. In this spirit, e.g., the 
famous Savvidy model , which favors a covariant constant magnetic field as ground state, 
has given rise to much speculation on the nature of the vacuum. 

Since a useful description of confinement and the ground state should also exhibit the 
limits of their formation, it is natural to perform a study at finite temperature where a 
transition to a deconfined phase is expected (as is observed on the lattice). An order 
parameter for the deconfinement phase transition in pure gauge theory is given by the 
Polyakov loop 0, E|, i.e., a Wilson line closing around the compactified Euclidean time 



1 



direction: 




L(x) = — trT exp \g \ dx k {x Ql x) , (1) 



where the period (3 — 1/T is identified with the inverse temperature of the ensemble in 
which the expectation value of L is evaluated. T denotes time ordering, N c the number 
of colors, and Ao is the time component of the gauge field. The negative logarithm of the 
Polyakov loop expectation value can be interpreted as the free energy of a single static color 
source living in the fundamental representation of the gauge group [|J. In this sense, an 
infinite free energy associated with confinement is indicated as (L) — > 0, whereas (L) ^ 
signals deconfinement. 

Moreover, (L) measures whether center symmetry, a discrete symmetry of Yang-Mills 
theory, is realized by the ensemble under consideration ||. Gauge transformations which 
differ at x = and x = (3 by a center element of the gauge group change L by a phase 
e 2mn/N c ^ n jaeger i eaV e the action invariant); this implies that a center-symmetric 
ground state automatically ensures (L) = 0, whereas deconfinement (L) ^ is related to 
the breaking of this symmetry. 

Therefore, the effective action governing the behavior of L is of utmost importance, 
because it determines the state of the theory at a given set of parameters, such as tem- 
perature, fields, etc. While this scenario has successfully been established in lattice for- 
mulations ||, several perturbative continuum investigations have led to various results. 
In the continuum, it is convenient to work with the "Polyakov gauge" , which rotates the 
zeroth component A of the gauge field into the Cartan subalgebra of SU(iV c ), A — > A 
(cf. Eq. (Q) below); furthermore, the condition 8qA = is imposed. Then, if the A 
ground state of the system is known, L can immediately be read off from Eq. (|1|), which 
suggests calculating the effective action of an time-independent A background field. 

Several one-loop calculations exist in the literature: considering a pure constant Aq 
background, Weiss || obtained an effective potential for the Polyakov loop preferring 
only the center-asymmetric ground states, i.e., the deconfined phase (see Eq. ([|D]) below). 
Combining the Savvidy model of a covariant constant magnetic background field with the 
Polyakov loop background Ao, Starinets, Vshivtsev and Zhukovskii as well as Meisinger 
and Ogilvie || were able to demonstrate the existence of a confining center symmetric 
minimum for (L) at low temperature with a transition to the broken, i.e., deconfined 
phase for increasing temperature. However, this model still suffers from the instabilities 
caused by the gluon spin coupling to the magnetic field f|, a problem also plaguing the 
Sawidy model ||, although an additional Aq background can in principle remove the 
problematic tachyonic mode in the gluon propagator for certain values of gAo and T. 



Perhaps the most promising approach was explored by Engelhardt and Reinhardt [ID|, 
who considered a spatially varying A field and evaluated the effective action for A in a 
gradient expansion to second order in the derivatives. The resulting action exhibits both 
phases, confinement and deconfinement, depending on the value of the temperature; in 
particular at low temperature, the spatial fluctuations of the Polyakov loop lower the action 
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when fluctuating around the center symmetric (confining) phase. The main drawback 
of this model is its nonrenormalizability. An explicit cutoff dependence remains; gauge 
and Lorentz invariance have been broken explicitly during the calculation. Nevertheless, 
the main lesson to be learned is that spatial variations of the Polyakov loop have to be 
taken into account while searching for an effective potential of the order parameter for the 
deconfinement phase transition. 

The present work is devoted to an investigation of the Polyakov loop potential partly 



in the spirit of flit)] ; however, the treatment of the quantum fluctuations, the calcula- 
tional techniques, and finally the results are quite different. In particular, we employ the 
background field method to keep track of the symmetries of the functional integral |IT] 



Unfortunately, the results are not as promising as those found in [|T(|, since the simple 



picture for the deconfinement phase transition is not visible in the most stringent version 
of the model. 

The paper is organized as follows: in Sec. |2|, we define the model, clarify our notations, 
and perform a first analysis of possible scenarios. Section [3] outlines the calculation of the 
effective action to one loop using the proper-time method, particularly emphasizing the 
subtleties of the present problem; we work in d > 4 dimensions with gauge group SU(iV c ). 
The implications of our results are discussed in Sec. [| for SU(2); therein it is pointed 
out that the main features of the model depend strongly on the treatment of the infrared 
sector. Section ||] briefly demonstrates the latter point by introducing an additional infrared 
scale "by hand" (gluon mass), which changes the properties of the model drastically, now 
exhibiting a confining phase. We finally comment on our findings in Sec. |[ 

One last word of caution: it is obvious from the very beginning that the one-loop 
approximation performed here is hardly appropriate for dealing with the strongly coupled 
gauge systems under consideration. In fact, the results presented below mostly represent an 
extreme extrapolation of perturbation theory to extraordinary large values of the coupling 
constant g without any reasonable justification. Nevertheless, besides being interesting in 
its own right, the model can serve as a starting point for more involved investigations. 
E.g., the renormalization group flow of the true effective action will coincide with the 
perturbative action at large momentum scales; hence, a detailed knowledge of the pertur- 
bative regime will be of use for checking nonperturbative solutions. Moreover, some of the 
technical results of the present calculation such as the form of the gluon propagator in a 
fluctuating A background will be expedient for other problems as well. 



2 The Model 

The essence of the model under consideration is determined by the choice of the background 
field, which is treated as a classical field subject to thermal and quantum fluctuations. At 
the very beginning, we confine ourselves to quasi-abelian background fields, pointing into 
a fixed direction n a in color space: 

k ll :=A a T a =:A li n a T a , n 2 = 1, (2) 
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where (T a ) bc = —if abc denote the hermitean generators of the gauge group SU(iV c ) in the 
adjoint representation. 

Now we are aiming at a derivative expansion in the time-like component^] of such 
an expansion is usually justified by demanding that the derivatives be smaller than the 
characteristic mass scale of the theory However, in the present case, there is no initial mass 
scale, since the fluctuating particles, gluons and ghosts, are massless. In fact, it turns out to 
be impossible to establish a unique derivative expansion for the (inverse) gluon propagator 
by a simple counting of derivatives; this is because a typical expansion generates terms 
~ ig.^4 | ; acting like a mass scale for the higher derivative terms. Therefore, we propose a 
different expansion scheme that is guided by the residual quasi-abelian gauge symmetry, 
which still holds for the background field. 

The model is further specified by assuming that there are no magnetic field components 
in the rest frame of the heat bath; the latter is characterized by its 4- velocity vector u^. 
Therefore, there are only two independent (quasi-abelian) gauge invariants: 

E 2 = -F^Fftv = F m u a FfupUp 



A„ 



- J dTA u (x" + tu^), A u :=A Pb u (i . (3) 



o 



Here we work in Euclidean finite-temperature space R d ~ l x S 1 , and F^ denotes the quasi- 
abelian field strength of the background field. In the heat-bath rest frame, we simply 
have u M = (1,0), so that A u = Aq. The quantity A u is invariant under quasi-abelian 



gauge transformations |L2| , since these transformations are restricted to be periodic in the 



compactified time direction. (For the complete gauge group, A u can be modified by a gauge 
transformation that differs at xo = and xq = (3 by a center element, e.g., A u — > — A u 
for SU(2) modulo Weyl transformations.) 

If we now perform a derivative expansion in the electric field E, we will obtain an 
effective action of the form Y = f(A u ) + g(A u ) E 2 + 0(E 4 ,Ed 2 E), (E = |E|, f,g to be 
determined) for reasons of gauge invariance. The indicated higher-order terms are at least 
of fourth order in dA and will be omitted in the following. 

Now, the crucial observation is that there exists a unique choice of gauge for the back- 
ground field that (i) satisfies the Polyakov gauge condition doA Q = in order to ensure 
the correspondence between Aq and L, and (ii) establishes a one-to-one correspondence 
between E and Aq, so that an expansion in E can be rated as a derivative expansion in Aq 

lr The spatially varying Ao field giving rise to an electric field appears to conflict with the assumption 
of thermal equilibrium, which is inherent in the Matsubara formalism used below; this is because electric 
fields tend to separate (fundamental) color charges, moving the system away from equilibrium. However, 
we adopt the viewpoint that the here-considered vacuum model characterizes only a few features of the true 
vacuum; the latter actually includes quark and magnetic gluon condensates (and higher cumulants) which 
altogether are in equilibrium. Beyond this, we expect the present approximation to hold for sufficiently 
weak electric fields, keeping the system close to equilibrium. Therefore, the expansion in the electric field 
performed below is consistent with (almost) thermal equilibrium. 
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(from now on, we work in the heat-bath rest frame where Aq = A u ): 

A (x) = a -(x- x')iEi Ei = -diA (x), (4) 

where a® and Ei are considered as constant, and x' is an arbitrary constant vector which 
can be set equal to zero. This gauge can be viewed as a combination of Polyakov and 
Schwinger-Fock gauge; the background field considered here lies exactly where the gauge 
conditions overlap. We remark that this is no longer true for higher derivative terms. The 
final task is to integrate out the thermal and quantum fluctuations in the background of 
the gauge field (Q) and expand to second order in E^. 

At this point, it is useful to introduce the dimensionless temperature-rescaled variable 

c:=||, CG [0,1]. (5) 

The compactness of c arises from the fact that Aq is a compact variable in finite-temperature 
Yang-Mills theories^. Then, the resulting effective action can be represented as a derivative 
expansion in c: 

r«Tff[e] = J d d x (V(c, n a ) + W(c, n a ) ftc^c) , (6) 

where the potential V and the weight function W also depend on the color space unit 
vector n a . Higher-order terms ~ (<9c) 4 are neglected. The superscript T in Eq. (^j) signals 
that the effective action strongly depends on the presence of a heat bath. Indeed, V as 
well as W vanish or reduce to simple constants as T — > 0; this is because the A field (or 
A u in Eq. (^)) ceases to be an invariant at T = 0, since it can be gauged away completely 
when the time direction is noncompact. 

Already at this stage, typical properties of the model become apparent. First, we 
observe that if W(c,n a ) > fluctuations of the Polyakov loop are suppressed; then, the 
ground state is solely determined by the minimum (or minima) of V(c,n a ), which we 
denote by cv,n,y. This ground state then is (not) confining if it corresponds to a center 
(a-)symmetric state, implying (L) = ((L) ^ 0). 

Fluctuations of the Polyakov loop can only be preferred if W(c,n a ) becomes negative 
for certain values of c and n a , which we denote by cw and n^. Whether or not these 
fluctuations lead to a confining phase again depends on the question of whether or not the 
minimum of W(c, n a ) corresponds to a center-symmetric state. Moreover, it depends on the 
question of whether these fluctuations are strong enough to compensate for the influence 
of V(c,n a ). Here we arrive at a main problem of the model: if W(cw,n^y) < 0, then the 
action is not bounded from below. In other words, arbitrarily strong fluctuations of c 
around cw will lower the action without any bound. Of course, it is reasonable to assume 



2 This can be inferred from a Hamiltonian quantization starting from the Weyl gauge Aq = and 
generating an Aq field by a time-dependent gauge transformation. This observation will furthermore 
become obvious when studying the background field dependence of the gluon propagator (cf. Eq. Jl^)). 
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(which we do in the following) that higher derivative terms (<9c) 4 or (d 2 c) 2 will establish 
such a lower bound, so that the strength of the fluctuations is dynamically controlled. 

Nevertheless, one drawback remains: we cannot make any statement about the nature 
of a possible phase transition. For this, we would have to know everything about the 
dynamical increase of dicdiC when W(cw, n^) becomes negative for certain values of tem- 
perature. Since this is beyond the capacities of our model, we shall always assume in the 
following that the system is dominated by the weight function W and thus by fluctuations 
of the Polyakov loop whenever W becomes negative. 

Let us finally perform a dimensional analysis of the model. For simplicity, let us start 
with d = 4. With regard to Eq. (^), the potential has mass dimension 4, while the weight 
function has mass dimension 2. Due to the compactness of A as reflected by Eq. (H), the 
only mass scale which is a priori present is given by the temperature T. Hence, if V scaled 
with T 4 and W with T 2 , say V(c, n a ) = T 4 v(c, n a ) and W(c, n a ) = T 2 w(c, n a ), where v, w 
are independent of T, then we would never encounter a phase transition in our model; 
this is because increasing or lowering the temperature could never turn W from positive 
to negative values or vice versa. 

At this stage, one may speculate that, since scale invariance is broken in Yang-Mills 
theories, the phenomenon of dimensional transmutation introduces another scale /i (e.g., 
the scale at which the renormalized coupling is defined). Then, the dimensionless function 
w can also depend on T / \l. However, this is far from self-evident, since the breaking of scale 
invariance is induced by UV effects. But the functions V and W in the effective action F^ s 
arise at finite temperature only and thus are a product of infrared physics. In particular, 
there are no UV divergences in the finite-temperature contributions to r e g- which require 
another scale during a regularization procedure. Hence, one is tempted to conclude that 
the naive scaling argument given above is correct. 

Nevertheless, the naive scaling breaks down, as we shall see in the next section; but 
this time, an additional scale is introduced by the properties of the theory in the infrared. 
As is well known, finite-temperature field theories can develop a more singular infrared 
behavior than their zero-temperature counterparts. Indeed, while the effective action at 
zero temperature and even the effective action for thermalized purely magnetic background 
fields do not suffer from infrared divergences, the case considered here involving thermalized 
electric fields exhibits such singularities, which must be handled carefully. The massless 
gluon does not provide for a natural infrared cutoff which could control the low-momentum 
behavior of the theory. 

To conclude, the d — 4 model is in principle capable of describing a phase transition, 
because the finite-temperature infrared divergences require an additional scale which in- 
troduces distinct high- and low-temperature domains. Of course, there are various ways 
to deal with the infrared singularities; and as we will demonstrate below, they can arise 
from different physical motivations, leading to different physical results. Two possibilities 
are proposed in the present work. In the first and more natural one, we regularize the 
infrared divergences in the same technical way as the ultraviolet ones, so that in toto there 
is only one more scale than in the classical theory, which we identify with the defining scale 
of the coupling constant. As a consequence and a consistency check, the running of the 
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coupling with temperature coincides with the running of the coupling with field strength 
or momenta - they are characterized by the same /3-function. In the second possibility, 
we study by way of example a regularization of the infrared divergences by an effective 
gluon mass m e g which we insert by hand, assuming that such an additional scale may be 
generated dynamically in the full theory. The latter version of the theory exhibits the 
desired properties of two phases separated by a deconfinement phase transition, while the 
former does not. 

At d > 4 the situation is somewhat different and simpler, since here the coupling 
constant g is dimensionful, so that two scales are present already at the classical level. 
Moreover, no additional scale will be introduced at the quantum level, because the theory 
is infrared finite for d > 4. 



3 Calculation of the Effective Action 

Starting from the standard formulation of Yang-Mills theories via the functional integral in 
Euclidean space with compactified time dimension, we employ the background field method 



11| to fix the gauge for the fluctuating gluon fields, but thereby maintain gauge invariance 
for the background field. We arrive at the one-loop approximation by neglecting cubic and 
quartic terms in the fluctuating fields. The remaining two integrals over the gluonic and 
ghost fluctuations are Gaussian and lead to functional determinants upon integration; the 
one-loop effective action depending on the background field then reads 

Tl s [A] = ±Tr xch In A™ [A]" 1 - Tr xc In A FP [A]' 1 , (7) 

where A ym [t4] _1 denotes the inverse gluon propagator, and A FP [v4] _1 the inverse ghost 
propagator, i.e., the Faddeev-Popov operator. The traces run over coordinate (x), color (c), 
and Lorentz (L) labels. Introducing the abbreviations D 2 := D^D^ and {DD)^ U := D fM D u , 
where the covariant derivative is defined by := — igA^, and suppressing the indices, 
the explicit representations of the propagators read 



AHA]" 1 



D' z -2igF+ I --1 } DD 



a 



Af[A]-i = —D 2 , (9) 

In the following, we will work in the Feynman gauge, a — 1, which simplifies the calcula- 
tions considerably^. For the evaluation of Eq. (^), the spectrum of the inverse propagators 
is required. In color space, diagonalization can be achieved by introducing the eigenvalues 
v\ of the matrix n a (T a ) bc , 1 = 1,..., N 2 — 1. The basic building block of the operators in 

3 For covariant constant background fields, the effective action is known to be independent of the gauge 
parameter pjJ. In the present approximation, the gauge field is covariant constant, so that we are 
allowed to set a = 1; this would no longer remain true if we were interested in the higher-derivative terms. 
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Eqs. ([8]) and (|9p is the covariant Laplacian, which upon insertion of the background field 
(|) yields (we set x' = 0) 

{-iD[A}) 2 = {-idif + (-iA)M) 2 + 2gu l E i (-iD [a ])x i + {gyi) 2 E l E j Xi x 3 , (10) 

where — i-Do[ a o] — — i<9o — 9 u i a 0i an d the roman indices run over the d— 1 spatial components. 
The operator is obviously of harmonic oscillator type and can be diagonalized by a rigid 
rotation of the spatial part of the coordinate system. A prominent eigenvector is given by 
the direction of the electric field E iy which we may choose to point along the 1 st direction 
of the new system. We finally obtain 

(—iD[A]) 2 = (-i<9 2 ) 2 + ■■■+ (-id^f + {gv l Ex 1 + (-iA))) 2 + (-idi) 2 , (11) 

where E = y/E^Ei. Up to now, we have achieved a partial diagonalization of the operators 
of Eq. (H) and @. While the Faddeev-Popov operator coincides with the Laplacian ([TT|), 
the inverse gluon propagator receives additional contributions from the gluon-spin coupling 
to the electric field ~ —igF^ v which can easily be diagonalized. Performing a Fourier 
transformation for the d — 2 unaffected components, — ic?2, • • • , — icta~i P2, ■ ■ ■ ,Pd-i, as 
well as for the time derivative, — id — > u n , uj n = 2nTn, n G Z (Matsubara frequencies), 
we may write the inverse gluon propagator in the form 

A™^]- 1 = pj + ■ ■ ■ + V \_ x + (gu l Ex 1 + (n )) 2 + (-i^) 2 + 2\gv l E, (12) 

where Uq = u> n — guia^. The number A labels the different eigenvalues in Lorentz space 
arising from the above-mentioned gluon spin coupling with A = 1, —1,0; here, A = 1, — 1 
appears only once, whereas A = occurs with multiplicity d — 2, corresponding to the 
spatial directions which are unaffected by the electric field. Incidentally, the Faddeev- 
Popov operator is identical to Eq. ( |T2"D with A = and multiplicity 1. Taking the prefactors 
and signs of the two traces in Eq. (0) into account, the Faddeev-Popov operator cancels 
exactly against two Lorentz eigenvalues of the spectrum of A YM L4] -1 with A = 0, removing 
the spurious gauge degrees of freedom, so that only the physical, transverse part of the 
inverse gluon propagator remains, 

Al^A]- 1 =p 2 2 + ---+ v\-x + (e, zi + (IIoh])) 2 + {-id x f + 2Xe h (13) 

where A = now occurs with multiplicity d — 4. For reasons of brevity, we introduced the 
short forms 

ei := \gUiE\, ai := |#v z a | (14) 

in Eq. (|I~3"D; the use of the moduli in Eq. ([14] ) is justified by the observation that, when 
tracing over a function of the inverse propagators, the result will not be sensitive to the 
signs of gviE and gvia . 

4 The compactness of A or ao becomes obvious here; e.g., for SU(2), where vi = —1, 0, 1, a shift of ao 
by an integer multiple of (2irT)/g can be compensated for by a shift of the Matsubara label n. 
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The remaining problem of diagonalizing the 0-1 subspace at first sight resembles the 
problem of finding the spectrum of a relativistic particle in a constant magnetic field. There, 
one finds the eigenvalues (Landau levels) by shifting the x% coordinate by x\ — ► X\ — ^"^^ 
in order to arrive at a perfect harmonic oscillator. Here, the situation is not so simple, 
because the ao field as well as the temperature dependence would drop out of the operator 
completely. In other words, such a shift is not in agreement with the periodic boundary 
conditions in time direction. 

Hence, the usual harmonic oscillator techniques arrive at their limits, and we have to 
find a different method that does not rely on the explicit knowledge of the spectra as 
is necessary for, e.g., ^-function methods. We choose Schwinger's proper-time technique, 
which provides for a more direct handling of the propagators. In terms of the transverse 
gluon propagator, the effective action reads in proper-time representation 

oo 

rl s [A] = ^Tr^lnAr^-^-itr^l^xle-^^-V) 

o 

oo 

- -l^J d i a t^ e '~'' (15) 



where Q denotes the spacetime volume of R d_1 x S 1 , s is the proper time, and the trace 
over the continuous part of the spectrum is taken in momentum space. The color trace 
runs over I, which labels the color space eigenvalues, whereas the Lorentz trace runs over A 
with its associated multiplicities. The function M is defined via the Fourier representation 
of the proper-time transition amplitude 



{x e 



sA™[A]~h n _sjh ip(x-x') -sM(p,X,l',») ( \Q) 

1 1 4^ (2vr) d ' K ' 



which can be determined by the differential equation 

1 = Al^A}- 1 A Y ± M [A}. (17) 
When evaluated, for example, in momentum space, Eq. (|P7D is solved by 

oo 

A Y ± M [A](p,X,l) = j dse- sM( P' x ' l ' s \ (18) 
o 

where M is given by 

„ w ,, x o o tanh2e;s. tanhe;s. , 9 
M{p,X,l;s) = Pl + ---+P 2 d ^ + — (Pi + qf + ^n-a* 2 

+ — lncosh2e ; s + 2ejA. (19) 

2t5 
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Here, q denotes some function of e\ and s which becomes irrelevant when shifting the pi 
integration in Eq. (|i"5|). Upon insertion of Eq. ([19]) into Eq. (|15|), the Gaussian momentum 
integration and the sum over A can easily be performed; the sum over Matsubara frequencies 
can be reorganized by a simple Poisson resummation,Q and we arrive at 



n 



-tr r 



(47r) d / 2 J S d l 2 




ds ( , d-2 
ei 4 sinners + — 



sinh eis 



00 ( n 2 
1 + 2^exp ( -—Qcoth 

n=l ^ 



< I* j cos— n 



(20) 



Here, we have separated the zero-temperature part, corresponding to the first line times 
the "1" of the second line, from the finite-temperature contributions, corresponding to the 
first line read together with the n sum. 



3.1 Effective Action at Zero Temperature 

Let us first study the temperature-independent part of the effective action Eq. (p0|) with 
particular emphasis on its renormalization: 

00 

rS-M = / + ■ < 21 > 



On the one hand, the proper-time integral is divergent at the upper bound, s — > 00, owing 
to the first term ~ sinhe/s. Since large values of s correspond to the infrared regime, this 
divergence is not related to the standard renormalization of bare parameters, which is a 
UV effect. In fact, this divergence is analogous to the Nielsen-Olesen unstable mode 
of the Sawidy vacuum^; one can give a meaning to this essential singularity by rotating 
the contour of the integral over the sinh term into the lower complex plane, —is — > s. The 
effective action then picks up an imaginary part that characterizes the instability of the 
constant electric background field considered here. 

On the other hand, the proper-time integral is also divergent at the lower bound, 
corresponding to the ultraviolet. The leading singularity is of the order s~ d ^ 2 , so that 
m subtractions are required for d = 2m or d = 2m + 1. The leading singularity which 
is field independent can easily be removed by demanding that r^f°[y4 = 0] = (first 
renormalization condition). The next-to-leading singularity proportional to ef ~ E 2 is 
removed by the second renormalization condition (dC e s/dE 2 )\ E _, = 1/2, where r eff = 
J £ e e, this ensures that the classical Lagrangian is recovered when all nonlinear interactions 

5 For technical details, see, e.g., Q ^4[ . 

6 At T = 0, the present situation involving an external electric field is identical to the magnetic Sawidy 
vacuum owing to the Euclidean 0(4) symmetry. 
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are switched off, and corresponds to a field-strength and charge renormalization 



L, n , = -rLrr> — -Z, a. 



f22) 



where .Er denotes the renormalized field, and Z 3 is the wave function renormalization 
constant. The latter can be read off from Eq. fl21|) by isolating the singularity ~ E 2 , 



_ 26 -d _ 2 
6(47r) d / 2 c9 



ds 



ad/2-1 



(23) 



where we have used an explicit cutoff A, employed tr c |z/;| 2 = \ u i\ 2 = N c , and have 

introduced the dimensionless coupling g 2 = g 2 n d ~ A with the aid of a reference scale fi (at 
which g is defined). To one-loop order, the (3 function can be read off from the coefficient 
of the UV divergence of : 



fa 



dfg 2 = (d - A)f - b d f 



(24) 



where 



60 



h d 



HA 

3 8tt 2 

(26 - d) N t 



3(d - 4) (47r) d / 2 



for d — i, 

for d > 4, 



(25) 



and t denotes the "renormalization group time" ln/x/A. Here we have rediscovered the 
standard well-known one-loop results, including the remarkable observation that the f3 
function for the dimensionful coupling g 2 = g 2 ji A ~ d vanishes precisely in the critical string 
dimension d = 26 [|T5 1 . 

Note that in 4 < d < 26, the j3 function develops a UV-stable fixed point: 



9* 



d-A 3(rf-4) 2 (47r) 



d/2 



h d 



26 -d 



(26) 



Of course, this fixed point lies in the perturbative domain (g 2 /47r 1) only for very large 
N c . 

As an alternative to these considerations of renormalization, the integral in Eq. (pi]) 
can be treated more directly with an appropriate regularization prescription. Let us briefly 
sketch a proper-time variant of dimensional regularization for later use in the case d = 4. 
Shifting the singularities at s — > by e and introducing a mass scale /i, Eq. ( |2"T| ) can be 
written as 



/*1T=0 
*-eff 



rtr c fx 



2c 



ds 

e 2-e 



e; sinh e«s + 



d- 2 



f 1 / 



s 2 ~ t sinhe^s 



(27) 
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These integrals can be evaluated [I6[, and the result for the one-loop contribution to the 
zero-temperature effective Lagrangian in d = 4 reads 



r^T=o |_ t 2 



f f 11 ej . . . 

— — in — + const. + imag. parts + Ofe) 

12e 12 ft 2 



The appearance of the simple pole in e implies a charge and field strength renormalization 
as outlined above. To be precise, in the background field formulation, the coupling runs 
with the scale set by the strength of the external field: g 2 = g 2 (gE/ /j, 2 ); this is analogous 
to the momentum dependence of the coupling in the standard formulation. Including the 
correctly (re-) normalized classical term, the total effective Lagrangian to one loop can then 
be written as 



C T £\gE) = £ cl + £1?=° = \b Q (gE) 2 In + imag. parts, re 2 = //e ^ \ (29) 

o ere 



where we have introduced the renormalization group invariant quantity re corresponding 
to the minimum of C^°(gE), and bo is given by the first line of Eq. (|2~5). 



Concerning the imaginary parts, the following comment should be made: within the 
Sawidy model, the imaginary parts indicate the instability of the constant-field vacuum 
configuration signaling the final failure of the model. In the present case, they are just an 
artefact of truncating the derivative expansion of the effective action at second order; this 
truncation is formally equivalent to the constant-field approximation. Upon an inclusion 
of non-constant terms which would affect only higher-derivative contributions, we expect 
the imaginary part to vanish; this is because the unstable modes are then cut off by the 
length scale of variation of the fields. 

For the regularization/renormalization program of £efr = ° in > 6, more subtractions 
than in d — 4 are needed; since these theories are nonrenormalizable in the common sense, 
these subtractions correspond to counter-terms of field operators of higher mass dimension. 
To be precise, d > 4 Yang-Mills theories can only be defined with a cutoff (with physical 
relevance); therefore, some cutoff procedure is implicitly understood. Nevertheless, the 
precise form of the cutoff procedure only affects higher-order operators which are of no 
relevance for the present model. Moreover, it is perfectly legitimate to study d > 4 quantum 
Yang-Mills theories in the sense of effective theories valid below a certain cutoff scale. 

3.2 Effective Action at Finite Temperature 

Equipped with these preliminaries, we now turn to the more interesting finite-temperature 
part of Eq. 

CO 

1T 1 f ds ( d — 2 \ ^-^ / n 2 \ a t 

CcS = "(47p^ tr V i^^^^^'+sinh^J 

(30) 



o 
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For s — > 0, the integral remains completely finite, since the coth in the exponent develops 
a 1/s pole; i.e., there are no UV divergences in the thermal contribution to £ e g, as is to 
be expected. 

At the opposite end, s — » oo, we again encounter the sinh divergence induced by the 
unstable mode. However, this is not the only infrared problem: an attempt at circumvent- 
ing this problem by a rotation of the s contour as in the zero-temperature case would lead 
to a disastrous behavior of the n sum due to the poles of the coth on the imaginary axis. 
In fact, it is the interplay between the proper-time integral and the n sum that produces 
further infrared divergences (at least for d=4). It is well known in the literature that par- 
ticles with Bose-Einstein statistics develop stronger infrared singularities at T 7^ than at 



zero temperature fll7| . Unfortunately, the status of these finite-temperature singularities 
is far from being settled, contrary to the T = case. 

In the present paper, we shall investigate two different methods. The consequences 
of an explicit mass-like cutoff are discussed in Sec. Here, we propose a more natural 
treatment by regularizing the thermal infrared divergences of Eq. (^) by the same method 
used to treat the UV divergences in Eq. ( ]2"7D in the T = case. Thereby, the same scale 
/j, which serves to define the value of the coupling constant is introduced. 

Taking these considerations into account, the Lagrangian is modified according to (sub- 
stitution fj, 2 s = u) 



00 



r>\T ■* . d I du { e l ■ r, e l , d ~ 2 e l/^ 



( 47r )d/2 J ud/2 - e y ^2 ^2 4 sinh £^ 





00 



00 / 2 2 \ 

^exp ^-IL^^coth^uJ cos|n. (31) 



x 

71=1 



In the context of our approximation in terms of derivatives of Aq, we need only the terms 
~ e° and ~ ef of Eq. (|31|). Expanding in e;//i 2 and performing the s integral, we arrive at 



7V r 2 -l 



- -™EE^^-mi (32 ) 



7T ' A — ' A — ' n 

1=1 n=l 



1T , i_ vV^-V v c -^> f^ 2+e 

eff| e 2 - Q d/2 2^ \jj2j 1^ n d I 4T 2 

1=1 VA ^ 7 n=l v 



X 



r(d/2-2-e) [(26 - d) - (d - 2){d - 4 - 2e)]T d , (33) 



For the term ~ e° in the first line, the e — > limit could safely be performed for d > 0; by 
construction, this term depends only on a\ ~ a ~ c (cf. Eq. @) and therefore corresponds 
to the potential V(c, n a ) as introduced in Eq. 

The term ~ ef in the second line contributes to the function W(c, n a ) (in addition 
to the classical term). It turns out that, for d > 4, the limit e —>■ can be performed 
immediately without running into an e pole. This means that, in these dimensions, the 
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thermally modified infrared behavior of the theory is under control. The order ef term of 
the one-loop effective action then reads 



C 1T 



T(d/2 



cS \ P 2 



g 67r d/2 



N£-l 

£ 

1=0 



2 • — COS if 11 



T 2 / * — ' w 

n=l 



d-4 



[(26- d) - (d-2)(d-4)]T d , d>4. (34) 



Obviously, the /x dependence has dropped out as a consequence of the well-behaved e — > 
limit. Nevertheless, there is a second scale besides the temperature, which is given by the 
dimensionful coupling constant g in d > 4. 

In d = 4, the situation is more involved, since Eq. (|3~3"D develops a simple pole in e for 
e — > 0. In order to isolate the pole and the terms of order e° which contain the physics, we 
first have to perform the n sum; this can be achieved with the aid of the polylogarithmic 
function (also Jongquieres function) 



Li(z,q) := 



71=1 



q" 

n z 



(35) 



and its analytical continuation for arbitrary real values of z []18 |. We finally find for Eq. fl33] 
in d = 4: 



V 2 -! 



/.IT I 



/=1 



11 11 T 2 11 11 -a; 

— - — In + — Li'(0, e 1 ^) + — Li'(0, e" 1 ^) 
12e 12 /i 2 6 V ' ' 6 V ; 



H C In 4 

6 12 12 



(36) 



where the prime at Li denotes the derivative with respect to the first argument, and C is 
Euler's constant C ~ 0.577216. Our first observation is that the e pole in this thermal 
contribution is identical to the one for the zero-temperature Lagrangian in Eq. (p8|). Since 
the latter is responsible for the usual charge and field strength renormalization leading to a 
field-strength-dependent coupling g 2 = g 2 (gE/ /x 2 ), the present e pole analogously suggests 
a running of the coupling with the scale set by the temperature: g 2 = g 2 (T 2 /fi 2 ). And 
because the residues of each pole are identical, the thermal running is governed by the same 
(3 function. This can be viewed as a consistency check of our treatment of the infrared 
singularities. 

Furthermore, the terms ~ e° depend on the ratio T 2 /ii 2 (even in the limit a\ — > 0). 
This implies that they cannot be normalized away as in the zero-temperature case, but 
lead to a thermal renormalization of the two-point function. 

This is in perfect analogy to QED, where an equivalent modification of the two-point 
function appears with the prefactor (= Yang-Mills j3 function) replaced by the QED (3 
function, and the role of /x is played by the natural scale of QED: the electron mass 

■0- 

In conclusion, it is the In 2^. term in Eq. (|36"D which leads to a breakdown of the naive 
scaling as outlined in Sec. ^and allows for a separation of high- and low- temperature 
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regimes. This could in principle facilitate a description of a phase transition within the 
d = 4 model. However, as we shall find in the next section, the model does not make use 
of this option. 



4 Analysis of the Effective Action 

In the following analysis of the previously derived effective action for arbitrary d and N c , 
for simplicity we confine ourselves to N c = 2, which provides for a convenient study of all 
the essential features of the model. Then, the color space eigenvalues V\ are simply given 
by 

14 = -1,0,1, for SU(2). (37) 

The results given above can be summarized in the effective Lagrangian (cf. Eqs. (||) 
and d)): 

£ e T ff [c] = V(c) + W(c) dicdiC, (38) 
where we have used the relations (cf. Eq. also (|)) 

c=^-, and <9 iC =^^, ce[0,ll. (39) 
2vrT' 2vrT ' L ' J V ; 

The convenient dimensionless quantity c is now considered as the dynamical variable of 
the effective theory; for SU(2), the center symmetric point is given by c = 1/2, since center 
symmetry relates c with 1 — c. If the vacuum state is characterized by c = 1/2, our model is 
confining, whereas a vacuum state different from c = 1/2 characterizes the deconfinement 
phase. 



4.1 Four Dimensions d = 4 

Beginning with the most relevant case of four spacetime dimensions, the potential can be 
read off from Eq. fl32|) . Performing the n sum leads to a Bernoulli polynomial, 

nc) = -^T 4 + ^TV(l-c) 2 , (40) 

in agreement with ||. While the first term is simply the free energy of N% — 1 = 3 free 
gluons, the second models the shape of the potential revealing a maximum at c = 1/2 
and minima at c = 0, 1 and thereby characterizing the deconfinement phase (see Fig. 
|l](a)). However, even if the potential had displayed a minimum at c = 1/2, it would have 
been of no use, since the potential by itself would remain confining for arbitrarily high 
temperatures. There would be no comparative scale separating two different phases. A 
Polyakov loop potential depending on c and T only can never model the deconfinement 
phase transition of Yang-Mills theories! 
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The weight function W(c) can be read off from Eq. (^) in combination with the classical 
Lagrangian £ d = E 2 /2 = ^^-dicdiC. 



W(c) = 2n 2 T 2 {^-b 



- In - + % + i- - In 2 + Li'(0, e 2 ™) + Li'(0, e' 2 ™) 

jJL I 11 



2% 2 T\ 



In X - I - 1 - ^ + In 2 - Li'(0, e 2 - ic ) - Li'(0, e" 2 ^) 



(41) 



where &o denotes the /3 function coefficient given in the first line of Eq. ( |25| ) (for N c = 
2). In the second line, we have expressed the running coupling and the scale fi by the 
renormalization group invariant k defined in Eq. (|2"D|), so that W(c) is itself renormalization 
group invariant! In fact, lowering the temperature can turn the weight function negative 
for any value of c so that fluctuations of the Polyakov loop are energetically preferred for 
low T. However, the confining value c = 1/2 always represents a local maximum of the 
weight function W(c), as is visible in Fig. |l|(b). For c — > 0, 1, the weight function diverges 
to — oo, but at c = 1, it jumps to its absolute maxima. Analytically, one finds 

W([c = 0, 1; c = 1/2]) = 2n 2 T% ([ln4vr; Iutt] - H - £ + l n 1\ . (42) 

To conclude, although our model indicates that fluctuations of the Polyakov loop become 
important at low temperatures, they do not fluctuate around the confining minimum, 
but energetically prefer a center asymmetric ground state for c. Hence, our model is not 
capable of finding a confinement phaseQ. Nevertheless, it should be stressed that the present 
treatment of the infrared modes is part of the definition of the model, although we have tried 
to formulate the present version as "universally" as possible. In fact, the regularization 
method considered here, which belongs to the standard class of regularization techniques, 
guarantees scheme-independent results. But it is also possible that the infrared modes 
are screened by a physical mechanism which involves another scale and thereby introduces 
"nonuniversal" information. Such a version of the model is discussed by way of example 
in Sec. ||. 

7 The discontinuous behavior of the weight function for c — > 0, 1 gives rise to speculations. Physically, 
such behavior is not acceptable (nor interpretable); rather, one may expect that some mechanism will lead 
to a wash-out of these singularities unveiling a smooth functional form of W(c) for c £ [0, 1] (although 
the origin of such a mechanism is still unclear to us). Probably, this will lead to a weight function of 
mexican-hat type with deconfining minima. However, with even more reservations, one might speculate 
upon the possibility of a smooth curve for W(c) which directly interpolates between the extremal values 
at c = 0, 1/2, 1 given in Eq. ( ff2| ) with a confining minimum at c = 1/2. Then, the model would exhibit 
a confining phase for small enough temperatures when W(c) becomes negative for c = 1/2. The reason 
for mentioning such vague speculations is to demonstrate how possible predictions could in principle arise 
from the model: following the reasoning of Sec. |[ the temperature of the phase transition is then given by 
W(c = 1/2)\t=t ct — 0. From Eq. (fl2|), we obtain: T CI /y/n~ 0.60. Identifying k with the string tension a 
(as it is the case in the leading-log model |2(J), our speculative estimate is in remarkably good agreement 
with the lattice value plf, T cr /^/a ~ 0.69. 
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Figure 1: (a) SU(2)-Polyakov loop potential V(c) in units of T d for d = 4, 7, 8, 9 (cf. 
Eqs. (fK]) and (f43f)). The d — 4, 8 potentials correspond to Bernoulli polynomials and 
Bg. (b) SU(2) weight function W(c) in units of k in d = 4 for different values of the 

temperature t := T/s/k = 0.2,0.6, 1 (cf. Eq. (f4l"l) ). The disconnected absolute maxima at 
W(c = 0, 1) are not depicted. 



4.2 Beyond Four Dimensions d > 4 

In spacetime dimensions larger than four, the situation simplifies owing to the absence 
of infrared problems. The Polyakov loop potential is again given by Eq. (^), which, for 
N c — 2, reads 



V{c) 



(d-2)r(d/2)c(d) 



7T 



d/2 



rpd 



d/2 



71 



(d 



. / » / \ COS 27TC7T, „ , 

2)r(d/2) ^ ^ T d , (43) 



n=l 



n" 



where ((d) denotes Riemann's ( function. Equation (^3|) is in perfect agreement with p2| 



where it is demonstrated that a representation of ^(c) in terms of Bernoulli polynomials of 

dth degree exists in d — 2, 4, 6, 8, We could as well choose a representation in terms of 

polylogarithmic functions which interpolate smoothly between the Bernoulli polynomials. 
In toto, the qualitative behavior of V(c) does not change significantly for different d: 
V(c = 1/2) is always a (deconfining) maximum (cf. Fig. |T](a)). 

The situation is different for the weight function W(c): in terms of the dimension- 
less coupling g 2 = fi d ~ A g 2 and polylogarithmic functions, the contributions from Eq. 
together with the classical Lagrangian can be represented as 



W(c) 



2 T 2 d _ 2 \ 1 T d -*T(d/2-2) 
/i 2 ^ I g 2 fj, d ~ 4 



48tt<V 2 



[(26-t/) - (d-2)(d-A)] 



Li(d - 4,e 



27ric\ 



Li(d - 4,e 



— 27ric\ 



(44) 



On the one hand, we again encounter the combination of polylogarithmic functions that 
interpolate between the Bernoulli polynomials of (d — 4)th degree for d = 6, 8, ... , es- 
sentially maintaining their typical shape. On the other hand, there is an important sign 
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Figure 2: (a) SU(2) weight function W(c) in units of // for d = 6, 7, 8, 10 and fixed T 
and g (cf. Eq. (f44l) ). Above d = d CT ~ 7.42, c = 1/2 represents the minimum of W(c). 
(b) The same weight function is now plotted for fixed g and d = 10 > d CT for various 
temperature values close to T cr . 



change owing to the factor (26— d) — (d— 2)(d— 4) at the "critical dimension" 

dcr= ^(5 + ^97) - 7.42. (45) 

For d < d CT , W(c) has a maximum at c = 1/2, implying that there is no confining phase in 
these dimensions. But for d > d cri the weight function exhibits an absolute minimum at 
the center symmetric value c = 1/2 (see Fig. 0(a)). As a consequence, W(c) can become 
negative at c = 1/2 for increasing temperature, as is depicted in Fig. 0(b). This is in 
agreement with the fact that the dimensionless coupling grows large in the /ug/i-momentum 
limit with a UV-stable fixed point given by Eq. (PB"|). Therefore, the model, somewhat 
counter-intuitively, describes a system with two different phases, a deconfined phase at low 
temperature and a confining strong-coupling phase at high temperature. In terms of the 
dimensionful coupling constant, the critical temperature where W(c = 1/2)\t=t ci — is 
given by 

q ^ = ^ tl d >d (46) 

9 cr r(d/2-2)C(d-4) (2 rf -5-l)[(d-2)(d-4)-(26-d)]' cr ' K } 

Because of the strong increase of the T and ( function in the denominator, the left-hand 
side rapidly falls off for increasing d. Typical values are g 2 T^~ 4 ~ 411.4,12.0,0.036 for 
d — 8, 16, 26. Therefore, the deconfined phase vanishes in the formal limit d — ► oo. 

Incidentally, it is interesting to observe that the discontinuities of the weight function 
W(c) for c — > 0, 1 vanish for d > 5; there, W(c) runs continuously to a finite extremal value 
for c — ► 0, 1. Between four and five dimensions, the discontinuity persists and W(c = 0, 1) 
increases for increasing d, finally approaching plus infinity at d — > 5~. 
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5 Additional Infrared Scales in d = 4 



The preceding section revealed that the d = 4 model required additional instructions on 
how to treat the singular infrared modes. Although we rate the procedure established 
above as the most general one of a "universal" character, we shall now suggest another 
method, involving an additional scale. In the following investigation, we exemplarily pick 
out one (physically motivated) possibility of regularizing the infrared modes, and study its 
consequences. 

Let us assume that Yang-Mills theory dynamically generates a scale in the infrared 
which can be reformulated in terms of an effective massQ m e g for the transverse fluctuating 
gluons[|. Although this scale may in itself depend on some parameters, we shall consider 
it to be constant within the limits of our investigation. 

Adding the effective mass term to the inverse transverse gluon propagator, e.g., in 
Eq. (0), it appears in a standard way in the proper-time representation of the effective 
action; for example, the integrand of the thermal one-loop contribution in Eq. (|30D is 

2 

multiplied by e~ m ^ s which damps away the infrared singularities. Upon the substitution 
u = m^gS, we obtain 

OO , s 

T _ _tt4 fdu ei ej 1 \ 

LcS ~ 4^ tlc J u* ml, [ Smh ml s U + 2 sinh ^-u ) 

00 / 2 2 \ 

El n*mi s ei e x \ a t 

exp -|2 — tt- coth — =-u cos — n. (47) 

Expanding for small eijm\ s in order to arrive at a consistent derivative expansion for A , 
we find to order e° and ef 

I °° T 2 ai 

^eff| = m cff tr c ^ — K 2 (m eS n/T) cos — n = V(c, n a , m cS ), (48) 

n=l 

II °° ai 
^eff| e 2 = "24^2 tr ° ^ ^2 K o( m esn/T) cos 

n=l 

+ 24^ tlc ^ ^ K i( m csn/T) cos ^n, (49) 

- n=l 

where we have employed a representation of the modified Bessel function K u (x) |I6| . Since 
we are interested in a possible formation of a confinement phase, let us study Eqs. ( [48]) 

8 This mass should not be associated with a thermal gluon mass; the latter represents a collective exci- 
tation of the thermal plasma and is a typical feature of the high-temperature domain, being proportional 
to T. By contrast, the effective mass considered here shall particularly affect the low-temperature modes 
and be approximately constant in T. 

9 In this way, gauge invariance with respect to the background field is maintained. 



19 



and ( fi^D in the low-temperature limit T <C m e fr- Then it is sufficient to use the asymptotic 
form of the Bessel functions for large argument, K v {x) — > a/vt /2xe~ x . 

Confining ourselves to the simplest case SU(2), we can deduce the form of the potential 
from Eq. (fi8|): 

V(c,m cS ) ~ -^Jt 4 (^y /2 e -m cB /T f cos27rc+ ^j j T « meff . (go) 

Again, we encounter a potential with a (deconfining) maximum at c = 1/2, so that the 
effective mass does not induce significant changes to the potential term. 

Including the contribution from the classical Lagrangian, the weight function can be 
deduced from Eq. ([4~9"D in the same limit: 



W(c,m eS ) = T 2 



2tt2 
9 2 




e -m cS /T cog 



(51) 



We first observe that, since the m e s dependent term is exponentially small for T <C m er T, 
a small coupling g 2 will always ensure that W(c, m e s) is positive so that Polyakov loop 
fluctuations are suppressed and the system is in the deconfined phase. Therefore, the 
model predicts that confinement requires a strong coupling. 

Indeed, if the coupling is (very) strong, we may neglect the first term in Eq. (0), and 
find that W(c,m e g) develops a minimum at c = 1/2, if0 

T < T cr , — = ^r, for<7 2 »l. (52) 
m eS 11 

The situation can be rephrased as follows: if T < T cr , c = 1/2 is the absolute minimum of 
W(c, m e fr). But W(c = 1/2, m e g) only becomes negative (thereby allowing for a confine- 
ment phase) if the coupling is sufficiently large, so that the first term of Eq. ([H]) can be 
neglected. 

Therefore, our main conclusion of the present section is that a different treatment of 
the infrared modes changes the behavior of the model significantly! Although the present 
version of the model exhibits the desired features, it requires more input and thus is less 
meaningful: we need to specify the value of m e a and the value of g 2 \ the latter is involved 
with another scale \i. 

Let us end this section with the comment that the introduction of a masslike infrared 
cutoff as employed in Eq. ( fTFD can also be used as an alternative regularization scheme for 
the infrared modes. This means that, giving up the meaning of m e g as a physical scale, 
but treating it as an arbitrary infrared cutoff scale for Eq. , we may remove it after the 
calculation by taking the limit m eS — > in Eqs. (|48|) and (f49"l). We exactly recover Eq. ( p2|) 
(for d = 4), and, after analytical continuation, also Eq. (^) with the association m e R ~ /i. 
The same procedure in d > 4 dimensions also leads to results identical with those in the 
preceding section. It is in this sense that the treatment of the infrared modes as performed 
in the preceding section can be rated "universal" . 

10 Taking the Bessel functions and the n sum more accurately into account, the actual value of T cr 
changes slightly: T cr /m e ff— 2/21. 
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6 Conclusions 



In the present work, we have established and analyzed a dynamical model for the order 
parameter of the deconfmement phase transition in Yang-Mills theories - the vacuum 
expectation value of the Polyakov loop operator. We have calculated the effective action 
for this order parameter to second order in a derivative expansion, and have treated the 
gluonic fluctuations in one-loop approximation. 

As a first conclusion, we observed that the "constant-Polyakov-loop" approximation, 
Aq =const., as considered in the literature, is in principle incapable of describing two differ- 
ent phases owing to the lack of an additional scale separating a high- and low-temperature 
phase in d = 4. This can also be inferred from the observation that the vacuum expec- 
tation value of the trace of the energy momentum tensor for a constant quasi-abelian 
background vanishes: 

(T%) = p g2 F^F^ = 0, for Al = n a A^ = const. (53) 

Therefore, a vacuum model of this type must necessarily preserve scale invariance even at 
finite temperature so that the theory must remain in a single phase. 

In the present model, scale-breaking is induced by fluctuations of the Polyakov loop 
which, in a particular choice of gauge, are associated with a nonvanishing electric field. 
The question of whether or not these fluctuations are energetically favored can in principle 
be answered by the dynamics of the model. In turns out that the deconfmement phase is 
the generic phase in the absence of fluctuations (this holds for all d > 4). Whether spatial 
Polyakov loop fluctuations drive the model into a confining phase depends on the form of 
the weight function W(c) of the kinetic term. 

In four spacetime dimensions, thermal infrared singularities complicate the investigation 
of the weight function and require additional specifications of how to deal with these 
singularities. Within a regularization-independent scheme that introduces no other scale 
than already present, the d — 4 model does not reveal a confinement phase; instead, 
fluctuations of the Polyakov loop even favor a deconfining vacuum state. 

By contrast, when regularizing the infrared by a physically effective cutoff comparable 
to a gluon mass for the transverse modes, a phase transition into a confining phase for low 
temperature becomes visible in the strong-coupling regime. 

Whether one of these scenarios is realized in Yang-Mills theory cannot, of course, be 
answered within a perturbative approach like that employed in the present paper. Not only 
does the enormous extrapolation of a one-loop calculation into the strong-coupling sector 
present a major problem, but, with regard to the infrared singularities, (even nonperturba- 
tively) integrating out the gluonic fluctuations at one fell swoop seems to be inappropriate. 
Instead, the integration over the fluctuations should be performed step by step in order to 
control a possible emergence of a dynamically generated mass scale. 

The one-loop model at least facilitates a concrete investigation of possible scenarios, 
and at most displays some features in a qualitatively correct manner. An appraisal of the 
different scenarios requires further arguments. The first scenario of Sect. ^J] without an 
effective mass can be preferred only from a theoretical viewpoint owing to its simplicity and 
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universality. Though the second scenario of Sect. || needs more input, the appearance of 
an additional infrared mass scale is common to almost all conjectured confining low-energy 
effective theories of Yang-Mills theory; therefore, a phenomenological viewpoint supports 
this scenario from the beginning, and so does the final result. Nevertheless, a reliable 
investigation of the infrared requires nonperturbative methods. 

Let us finally comment on the differences of our results to Ref. |K| which inspired 
the model considered in the present work; although the representations of the effective 
action in the form of Eq. are congruent, the meaning of the results is quite different 
in comparison: in [ lOfl , the fluctuations of the Aq field have not been taken into account, 



implying that the resulting "effective action" remains a complete quantum theory of the 
A Q field. The Aq ground state is then approximated by the effective potential which is 
obtained by transforming the kinetic term to standard canonical form. By contrast, we 
integrated over all quantum fluctuations of the A^ field in the present work; therefore, the 
resulting effective action is the generating functional of the 1PI diagrams and governs the 
dynamics of the background fields in the sense of classical field theory. To conclude, it is 
not astonishing that the explicit results of in particular for the weight function W (c) 
do not agree with ours because they have a different origin and a different meaning.^ 

In spacetimes with more than four dimensions, the situation simplifies considerably: 
on the one hand, the model is infrared finite, thereby producing unambiguous results; 
on the other hand, there already exists another dimensionful scale given by the coupling 
constant. We discovered a phase transition from the (generic) deconfining to a confining 
phase for increasing temperature for d > d cr ~ 7.42. This is consistent with the fact that 
the dimensionless coupling constant grows for increasing energies, reaching an UV-stable 
fixed point. Beyond perturbation theory, the latter statement has also been confirmed in 
the nonperturbative framework of exact renormalization group flow equations 
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